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MOTIVATION 


• Unsteady CFD flow calculations are computationally expensive when 
compared to steady flow calculations 

• Conflicting interests: We want adequate spatial & temporal accuracy 
but we don’t want to pay the price (Excessive CPU time) 

• The computational mesh drives the cost of CFD calculations and 
should be optimized for each flow condition. This suggests that 
solution algorithms should be closely tied with grid generation 

• How do we optimize the mesh? Distribute the numerical error evenly 
throughout the mesh 

• Use adaptive meshing to evenly distribute the spatial discretization 
errors 

- locally enrich in regions of relatively large errors 

- locally coarsen in regions of relatively small errors 
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OVERVIEW OF 2D MESH ENRICHMENT STRATEGIES 


• Type-4 enrichment element 



• Type-2 enrichment element 
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OVERVIEW OF 3D MESH ENRICHMENT STRATEGIES 







DESCRIPTION OF 2D & 3D UPWIND-TYPE EULER ALGORITHM OF BATINA 

• Finite-volume spatial discretization on unstructured-grids 

- triangles in 2D 

- tetrahedra in 3D 

• flux vector splitting of van Leer 

• Flux limiting to suppress oscillations near shock waves 

• Time integration may be either explicit Runge-Kutta scheme or implicit 
Gauss-Seidel relaxation scheme 

• Implicit scheme allows very large CFL numbers for rapid convergence 
to steady state 

• Choose time step for unsteady calculations based on physics of 
problem rather than numerical stability 






INSTANTANEOUS MESH AND DENSITY CONTOUR LINES FOR THE 
SHOCK DIFFRACTION PROBLEM 


• t = t 2 



Rausch, 1992 X * W ^ 








INSTANTANEOUS MESH AND DENSITY CONTOUR LINES FOR THE 
SHOCK DIFFRACTION PROBLEM 

• M s = 2.81 

• Ap = 0.2 


• t = t 3 • t = t 4 



INSTANTANEOUS MESH AND DENSITY CONTOUR LINES FOR THE 
SHOCK DIFFRACTION PROBLEM 


• M s = 2.81 

• Ap = 0.2 


• t = t 5 • t = t 6 








COMPARISON OF SHOCK TRIPLE POINT LOCATIONS WITH 
EXPERIMENTAL DATA 


• Experimental data by Bryson and Gross, Journal of Fluid Mechanics, 
vol. 10, pp. 1-16, 1961 
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PARTIAL VIEW OF THE SURFACE MESHE FOR THE SYMMETRY PLANE 
AND THE ONERA M6 WING 


• Total mesh has 46,516 tetrahedra and 8,824 nodes 
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COMPARISON OF UPPER SURFACE DENSITY CONTOUR LINES 
FOR THE ONERA M6 WING 

• M„= 0.84, a 0 = 3.06° 


• Ap = 0.025 


• Original mesh • 1 level 

• 2 levels 
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COMPARISON OF COEFFICIENT OF SURFACE PRESSURE 
FOR THE ONERA M6 WING 


• Original mesh • 1 level • 2 levels 
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ILLUSTRATION OF THE SHOCK-TUBE PROBLEM 


(1) High pressure 


(4) Low pressure 
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SURFACE MESH FOR THE SHOCK-TUBE PROBLEM 


• Total mesh contains 562 nodes and 1 ,800 tetrahedra 
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COMPARISON OF THE VARIATION OF DENSITY, VELOCITY, AND 
PRESSURE THROUGHOUT THE SHOCK-TUBE 


• Solution at time t = 0.1 

• Original mesh • Adapted mesh 



COMPARISON OF THE VARIATION OF DENSITY, VELOCITY, AND 
PRESSURE THROUGHOUT THE SHOCK-TUBE 


• Solution at time t = 0.3 

• Original mesh • Adapted mesh 









SUMMARY 


Final solution adapted mesh depends on the original mesh 

- adapted mesh cannot be coarser than the original mesh 

Enrichment/Coarsening procedures are robust for isotropic cells; 
however, enrichment of high aspect ratio cells may fail near boundary 
surfaces with relatively large curvature 

Enrichment indicator worked well for the cases shown, but in general 
requires user supervision for a more efficient solution 



